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Abstract 


Steam reforming of methanol is investigated numerically considering both heat and mass transfer of the species in a packed bed microreactor. 
The numerical results are shown to be in good agreement with experimental data [M.T. Lee, R. Greif, C.P. Grigoropoulos, H.G. Park, F.K. Hsu, 
J. Power Sources Transport in, 166 (2007) 194-201] with a BASF F3-01(CuO/ZnO/A1,03) catalyst. A correlation for the conversion efficiency 
of methanol has been obtained as a function of the operating temperature and a dimensionless time parameter which represents the ratio of the 
characteristic time of the methanol flow to the time for chemical reaction. The results show that for the constant wall temperature condition the 
steam reforming process of methanol results in a nearly uniform temperature throughout the microreactor over the range of operating conditions. 


Published by Elsevier B.V. 
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1. Introduction 


There is great interest in producing hydrogen for small 
fuel cell applications to provide high volumetric and gravi- 
metric energy density for miniaturized portable electronic 
systems [1]. Of special interest is the need to develop effi- 
cient microreformers to convert hydrocarbon fuel into hydrogen 
which can be delivered to a proton exchange membrane fuel 
cell to produce electricity. Many investigations indicate that 
methanol is an attractive fuel since it has a high hydrogen- 
to-carbon ratio relative to gasoline and lower inter-carbon 
bonds; hence it can be reformed efficiently at moderate tem- 
peratures (200-300 °C) [2]. The methanol reformer operating 
at moderate temperature also produces low CO concentra- 
tion which is important because of CO contamination of the 
anode. 

Our objective is to examine the conversion efficiency for the 
production of hydrogen and CO from methanol in a catalyst 
reforming process. Steam reforming of methanol is achieved 
by the chemical reaction over a catalyst packed bed at mod- 
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erate temperature. The steam reforming process of methanol 
is characterized by the operating temperature, velocity of the 
steam—methanol mixture reactant, the mass of the catalyst, 
etc. 

Many studies have been carried out to provide insight 
into the thermal performance and the heat transfer in a 
steam—methanol reforming reactor. A number of kinetics mod- 
els for steam—methanol reforming have been reported in 
the literature [3-8]. Amphlett et al. [3] suggested a semi- 
empirical kinetics model for the steam—methanol reforming 
process over the CuO/ZnO/A1203 catalyst, which includes the 
mathematical model of the reaction rates for the plug flow 
reactor. Park et al. [9] carried out a one-dimensional anal- 
ysis for the mass transport in a steam—methanol reformer 
that employs Amphlett’s reaction kinetics model. Karim et 
al. [10] carried a two-dimensional pseudo-homogeneous anal- 
ysis for packed bed reactors. In this model, they took the 
diffusion coefficient of the gas mixture in the reaction pro- 
cess as constant and focused on the effect of heat transfer 
on the measured methanol conversion. Park et al. [11] car- 
ried out a quasi three-dimensional analysis of the reacting flow 
in a steam—methanol reformer, coupled to three-dimensional 
heat transfer in a silicon wafer. Yuan et al. [12] simulated in 
a three-dimensional study the heat and mass transport in a 
methane-steam reforming duct considering the catalytic chem- 
ical reactions. 
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Nomenclature 


Porous surface area per unit volume of pellet 
[m? m~*] 

Cross-sectional area of bed [m7] 
Pre-exponential term in the Arrhenius expression 
for k 

Concentration of methanol [mol m~?] 

Specific heat of species i [Jkg~! K] 
Modification factors of decomposition and 
reforming reaction 

Diameter of catalyst pellet [m] 

Characteristic diffusion length [m] 

Effective diffusivity [m* s~!] 

Ordinary diffusivity [m° s~!] 

Knudsen diffusivity [m2 s7!] 

Activation energy in the Arrhenius expression for 
k [Jmol~!] 

Feed rate of liquid mixture [m3 s7!] 

Thermal conductivity of species i [W m7! K] 
Volumetric reaction constant for decomposition 
{mol kg s7 1] 

Modified reaction constant for decomposition, 
(1 — £)psk{} [mol m~? s] 

Surface reaction constant [m s7 !] 

Volumetric reaction constant for reforming 
[m3 kg! s] 

Modified reaction constant for reforming at 
T=220°C, (1 — e)psky [s7!] 

Enthalpy of reaction [J mol7!] 

Length of bed [m] 

Lewis number, Le = pgCp.gD m/kg 

Mass fraction of species i 

Mass of catalyst [kg] 

Volumetric mass source [kg m`? s] 

Cross sectional area-weighted average of CO 
mass concentration [ppm] 

Molecular weight of species i [kg mol~!] 
Pressure [Pa] 

Volumetric heat source [W m7?] 

Radial and axial Cartesian coordinates [m] 
Volumetric chemical reaction rate of species i 
[mol m7? s] 

Pore size of pellet [m] 

Universal gas constant [J mol`! K] 

Reynolds number 

Surface area [m?] 

Schmidt number 

Sherwood number 

Molar ratio of steam to methanol 

Pellet surface [m?] 

Temperature [K, °C] 

Velocity of mixture [m s71] 

Volume [m?] 

Mole fraction of species i 


Xi Cross sectional area-weighted average of mole 
fraction of species i 


Greek symbols 
Porosity 
Conversion efficiency 
Dynamic viscosity [Pas] 
Viscosity [m2 s7!] 
Density [m? kg !] 

52 Lennard Jones Parameters 
Tortuosity factor 


aqad IES % 


Subscripts 

Inlet, reference 
Bed 
Decomposition 
Gas mixture 
Species 

Liquid 
Average 

Pellet 
Reforming 
Solid 


opURrC TR gro 


In this work, we have analyzed the methanol—-steam reformer 
considering the two-dimensional variation of the transport and 
thermal properties of the gas mixture associated with the 
composition of each species throughout the reformer. The 
effect of wall temperature on the conversion efficiency and 
CO production from methanol in the reformer has been inves- 
tigated with the feed rate of the water-methanol mixture 
and the mass of the catalyst as parameters. For the com- 
mercial BASF F3-01(CuO/ZnO/A1203) catalyst, a correlation 
for the conversion efficiency of methanol has been presented 
as a function of the wall temperature of the reformer and 
a dimensionless time parameter that is related to the flow 
rate of the steam—water mixture fed to the reactor, the cat- 
alyst mass and the reaction rate. In addition, the numerical 
results for the conversion efficiency of methanol and the amount 
of CO produced in the reformer have been compared with 
experimental data [13] and with a one-dimensional analytical 
solution [9]. 


2. Analysis of steam methanol reforming process 


2.1. Chemical kinetics of catalyst steam methanol 
reforming 


Methanol can be reformed by two overall reactions in a 
reformer filled with the catalyst CuO/ZnO/Al203 as described 
by Amphlett et al. [3]: 

CH30H + H20<8+CO» + 3H2 a) 


CHOH < CO + 2H (2) 
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where the subscript ‘R’ refers to steam reforming, and ‘D’ 
refers to decomposition. Methanol is primarily reformed via 
reaction (1) called reforming. Some portion of the methanol 
also decomposes directly to CO via reaction (2). In addition, 
the water—gas shift reaction adjusts the composition of the 
product gas: 


CO + H20 <| CO2+ H2 (3) 


The water—gas shift reaction [reaction (3)] can be neglected 
without a substantial loss in accuracy of methanol concentra- 
tions as described by Amphlett et al.[3]. For a packed bed, the 
chemical reaction rates can be expressed in terms of reaction 
rate constants and the concentration of methanol according to 


rR =(1— 8) pskgci, rh = (1 — 8) pskp (4) 


where p; is the solid density of the catalyst, and ky’ and kğ 
are the rate constants for reforming and decomposition, respec- 
tively. The reaction rate constants depend on the properties of 
the catalyst surface and on the temperature. These quantities can 
be related to the rate constants of Amphlett et al. [3] which are 
derived from an Arrhenius relation as follows: 


ki’ = Cr(Ap + Br In SMR) (28/87) (5) 
kË = CpApe™ 20/R?) (6) 


where SMR is the molar ratio of steam to methanol; Apr, Br, and 
Ap are Amphlett’s constants ([3]) for reforming and decompo- 
sition reactions, respectively. Cr and Cp are correction factors 
for reforming and decomposition reactions, respectively, and 
can be determined empirically from the activity and effective- 
ness of catalyst. From the experimental data [13] for BASF 
F3-01(CuO/ZnO/Al203) catalyst, the correction factors are 
determined to have the values of 5.5 in the reforming reaction and 
3.5 in the decomposition reaction, respectively. The generation 
rates of the species can be expressed as 


rt =r — rg (7a) 
r3 = rp (7b) 
r3' = rR (7c) 
ry =p (7d) 
r! = 3rg + 2rp. (7e) 


where the subscripts are 1 = CH3OH, 2 = H20, 3 = CO7, 4 = CO, 
and 5 = H3. 


2.2. Heterogeneous catalytic reaction rate 


The methanol reforming reaction occurs on the catalyst par- 
ticles in a packed bed during the reforming of a steam—methanol 
mixture. It also includes the heat and mass transfer in the cat- 
alyst porous medium of the packed bed. The mass transport 
considerations control the rate of reaction, while heat transfer 
considerations control the rate at which the heat of reaction can 
be removed, and hence the temperature of the reaction. Our 
analysis for reaction rate is undertaken by considering the mass 


diffusion of methanol from the free stream to the surface of the 
catalyst pellet, and into the pores of the pellet. The heteroge- 
neous reaction rate within the packed bed of catalyst rg’ can be 


expressed by [14,15] 
(8) 


where Sp is the surface of pellet within the packed bed of cata- 
lyst, Vp is the volume of the pellet, and rg is the molar flux of 
methanol at the surface of pellet. The methanol-steam mixture 
diffuses into the pores of the pellet and encounters two resis- 
tances; one is a convective resistance between the free stream 
and pellet surface, and the other is a diffusive resistance within 
the pellet pore. The molar flux rg depends on both the con- 
centration of the methanol at the surface of pellet and the two 
resistances [16] 


cl 


id 
R= 


Dr (9) 
ShDim 


1 
(Vpap/Sp) pkg 


where cı is the concentration of methanol in the free stream 
between the pellets, ap is the porous surface area per unit vol- 
ume of pellet, Dp is the characteristic diffusion length in the 
packed bed defined by Dn = £€dp/(1 — £), np is the effectiveness 
of catalyst pellet, and kg is the first order kinetic rate con- 
stant for the reaction. The kinetic rate constant kg in Eq. (9) 
is defined per unit the surface area of pellet and related to 
the reaction rate constant kg in Eq. (5), kg = kg Ps/ap. The 
effectiveness of the catalyst pellet is the ratio of the actual 
reaction rate divided by that for a pellet with an infinite dif- 
fusion coefficient and can be taken for a spherical pellet to 
be [16] 


1 1 l v, [e 
tp = A= 2, 2R (10) 
A \tanh3A 3A Sp \ Di,eff 


where Dy eff is the effective diffusivity within the catalyst pellet 
and includes ordinary diffusion (Dım) and Knudsen diffusion 
(Dx1) [16]: 


L tee a 1 
= + (11) 
Dietf €p \Dim Dki 
1-x 
lm = a A.) 
YD peo(xi/Dii) 
T? Gh + in) 
Dii = 1.86 x 107’ +—+—___ 12 
2 SRT 
Dy = 3 pore aM, (13) 


The correlation for a dry packed bed for the Sherwood number 
(Sh), defined as a dimensionless mass transfer conductance in 
Eq. (9), depends on Reynolds number (Re) and Schmidt number 
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(Sc) as follows[17]: 


Sh = (0.5 Re0.5 + 0.2 Re2/3)Sc!, 


UD, 
= ial aed Sc= a (14) 
Mg Dim 


Re 


where the density of the gas mixture pg is the mass-weighted 
average of the density of species assuming an ideal gas mixture 
and can be expressed as follows: 


P 


P 
P= eM = Rp (15) 


The viscosity of the gas mixture 4g in Eq. (14) can be calcu- 
lated from Wilke’s mixture rule [18] as follows: 


5 
Xi fi 
a us 
ial jel iXij 


y, — + e/a My AP 
: V8 [1 + (Mi/M p"? 


(17) 


2.3. Conservation equations 


The methanol-steam reformer configuration is axisymmetric 
as shown in Fig. 1, and includes the isotropic and continuously 
distributed catalyst. The flow of the steam—methanol mixture 
is assumed to be steady and radially uniform through the inert 
catalyst packed bed. The mixture consists of five species, taken 
to be ideal gases, and reacts heterogeneously with the catalyst. 
It is also assumed that the thermal state of the catalyst and 
the gas mixture is locally in equilibrium and the Lewis num- 
bers of all species are equal to unity, Le=pgcp oD m/kg = 1. 
The Cartesian cylindrical coordinate system is used with the 
z-axis coincident with the apex of the flow region as shown 
in Fig. 1. The mass conservation of the gas mixture can be 
expressed as 


d 

—(epgU) = 0, 18 
az! PgU) (18) 
where U is the mean velocity of the gas mixture, and € is the 


catalyst porosity given by e= 1 — (V/V). The equation for the 
conservation of energy reduces to [19] 


aria aT\ a aT 
EpeCp,gU— = (rn )+ (rn ) +d (19) 


əz rðr ar az az 


Ly 


Fig. 1. Schematic of the modeling domain in the reactor. 


where km is a volume-weighted average conductivity of the cat- 
alyst material and gas mixture, and q/" is the volumetric heat 
source of the porous medium generated from the heterogeneous 


reaction between the gas mixture and catalyst: 


km = (1 — £) ks + ekg. (20) 
di, = —AHprR — AHprp. (21) 


The specific heat of the gas mixture cp» is obtained from 
the mass-weighted average for each species, Cp,g = Sm iCp,i 
[16]. The thermal conductivity of the gas mixture kg is calculated 
in the same manner as Hg (cf. Eqs. (16) and (17)) from Wilke’s 
mixture rule [18]. The conservation of mass equation for species 
i reduces to [19] 


Oe 12 rop ZEN a 2 op E N q m”. 
Pg əz rər Pg ™ or az Ps ™ az ei 


(22) 


where m; is the mass fraction of gas component i defined 
by m;=;/pg, Dm is the mean diffusivity of the gas mixture 
in the porous media, and my. ; is the mass generation of gas 
component i: 


Dm = €D1m. (23) 
mej = Mir,’ (24) 


A constant tube wall temperature, T= Ty, is specified. The 
temperature of the gas mixture at the inlet is taken to be the same 
as that of the wall. The mass transfer for each of the species is 
impermeable through the tube wall, 0m;/dr|,—a,/2 = 0. The liq- 
uid mixture of water and methanol is considered to have been 
preheated and then flows as a gas mixture into the packed cata- 
lyst at z=0. The molar ratio of water to methanol is prescribed 
and the velocity of the gas mixture is uniform. The physical 
parameters and basic operating conditions are listed in Table 1 
and the other physical gas properties of each species are referred 
to the existing literatures [16,20—22]. 


2.4. Numerical procedure 


The governing equations were solved numerically using the 
method of Karki and Patankar [23]. The governing equations 


Table 1 

Geometric parameters and physical properties 

Diameter of catalyst bed (dp) 1.0 x 10-3 m 
Axial length of catalyst bed (Lp) 10 x 10-7 m 
Diameter of pellet (dp) 150 x 107 m 
Porous surface area of pellet (ap) 1.92 x 108 m? m~? 
Pore size of pellet (7pore) 1.0 x 1077 m 


Porosity of catalyst bed (e) 0.35 

Density of catalyst (ps) 1300 kg m~? 

Thermal conductivity of catalyst (ks) 20 W m~! K (Touloukian [22]) 
Void—tortuosity ratio factor (€p/Tp) 0.2 (Mills [16]) 

Lennard Jones parameter (01;) (Mills [16] A.26) 

Lennard Jones parameter (2p) (Mills [16] A.28) 

Molar ratio of steam to methanol (SMR) 1.1 
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Fig. 2. Molar fraction variations of the components for the feed water-methanol flow rate of 10 pl min™!: (A) molar fraction variations of species along the axial 
direction of the reactor and (B) molar fraction variations of species with respect to Ty leaving the reactor. 


are transformed into the curvilinear cylindrical coordinate sys- 
tem r=r(&, n) and z=z(&,n) and the equations are discretized 
using the central difference scheme [24]. The resulting alge- 
braic equations were solved using a nonuniform grid system 
with 71 nodes in the é direction and 45 nodes in the 7 direction. 
In the region near the surface and in the vicinity of the tube 
wall, the grids were densely distributed. Iterations were contin- 
ued until changes in the mass fractions were less than 0.1%. It 
was found that the results for the temperatures and mass frac- 
tions are differed by less than 0.01% from a grid system of 
90 x 142 nodes. 


3. Results and discussion 


In this work, results are obtained for a methanol 
steam reformer for a circular tube packed with the 
BASF F3-01(CuO/ZnO/A1203) catalyst, and investigated over 
the following range of conditions: 200°C < Ty < 240°C, 
Smg < mg, < 25mg, 5 pl/min < Fjo < 25 ul min—!, 0.5mm < 
dy < 4mm, 2mm < Lp < 20mm. 

Fig. 2 shows the molar fraction variations of the components 
produced for the feed water—methanol flow rate of 10 pl min7!. 
As shown in Fig. 2(A), when the steam—methanol mixture 
decomposes at the constant wall temperature, Tw =220°C, the 
main components produced are hydrogen and CO2, and only a 
small amount of CO is produced, which is not discernible in the 
figure. As the steam—methanol mixture is gradually reformed 
throughout the reactor, the molar fractions of hydrogen and 
COz increase along the axial length of the reactor. Fig. 2(B) 
shows nearly linear variations of the mole fractions of each 
species leaving the reactor as a function of the wall tempera- 
ture. The increase of Hz produced, as shown in Fig. 2(B), in 
the reformer is attributed to the activation of the chemical reac- 
tion of methanol with respect to an increase of the operating 
temperature. 


The conversion efficiency of methanol as a function of the 
wall temperature is shown in Fig. 3 for three values of the cat- 
alyst mass for a feed water-methanol flow rate of 10 wl min7!. 
The effect of increasing wall temperature is to almost linearly 
increase the conversion efficiency of methanol in the reformer. 
The methanol conversion rate increases with increasing catalyst 
mass. 

The influence of the catalyst mass on the flow rate of the 
products leaving the reactor is shown in Fig. 4 for the mixture 
feed flow rate of 10 ul min~!. The molar flow rates of Hy and 
COz increase with increasing catalyst mass. The flow rate of CO 
also increases with the catalyst mass, but is not discernible in 
the figure. 


a 
a 
a 
„e 
a 
_— 
— 
—_ — q 
210 220 230 240 
T,(°C) 


Fig. 3. Conversion efficiency of methanol for three cases of catalyst mass for 


the feed water-methanol flow rate of 10 pl min™!. 
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Fig. 4. Influence of the catalyst mass on the flow rate of species leaving the 


reactor for a feed water—methanol flow rate of 10 wl min™!. 


The conversion efficiency of methanol for three values of the 
feed flow rate of the mixture is shown in Fig. 5 as a function 
of the wall temperature. The conversion efficiency of methanol 
decreases with increasing feed flow rate. This is because the 
mixture stays in the catalyst packed bed for a shorter time with 
increasing flow rate, which results in less reaction with the cat- 
alyst. 

Peppley et al. [4] refer to a fractional methanol conversion 
as a function of a parameter related to the ratio of the catalyst 
mass to the feed flow rate of the mixture, ms/F\,9. This parameter 
can be modified into a dimensionless parameter, msky o/ Fi,0Ps- 


1.0 T T T T T T T 


F,o (#l/min) 


1 l 1 1 1 1 1 
°-So0 210 220 230 240 


Fig. 5. Conversion efficiency of methanol for three cases of feed flow rate for 
ms = 16mg. 
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Fig. 6. Results for the conversion efficiency of methanol as a function of the 
dimensionless parameter. 


Considering the relations ms = psAbLp and Fi o= UA», this 
parameter is related to the parameter, Lpk’”R,o/ U(cf. Eq. (27) 
discussed below, where kg o = (1 — £) pskp’ o). This dimension- 
less parameter is related to the ratio of the characteristic time 
of the methanol flow, m;/Fi,ops to the time for the methanol 
chemical reaction in the reformer, 1/ kg o- Results for the 
conversion efficiency of methanol are shown for four wall tem- 
peratures in Fig. 6 in terms of the dimensionless time parameter, 


msky o//i,0Ps- The conversion efficiency of methanol increases 


for increasing values of mskg o/F1,0ps. Based on the numerical 
results, the conversion efficiency of methanol is correlated as 
follows: 


n = —1.222 + 0.00061T (K) 


Ala 


mskp o (25) 
Fi, 0Ps 


— (1.087 — 0.0033T (K)) In 


This correlation agrees with the calculated conversion effi- 
ciency with an error less than 0.1% over the range of operating 
conditions. 

The numerical results for the conversion efficiency of 
methanol are now compared with both experimental data [13] 
and the results from a one-dimensional analysis. Considering 
the molar ratio of steam to methanol, SMR, the one-dimensional 
analytical solution for the conversion rate of methanol is given 
by [9,13]: 


me ŤCH30H,0 — XCH30H _ E(SMR + 3) (26) 
7 XCH30H,0 2§€+SMR+1 
where é is 
L 
2&+(SMR+ 3)In(l—&) = 9 (27) 


U/[K"R(SMR + 1)] 


J.-S. Suh et al. / Journal of Power Sources 173 (2007) 458-466 


1.0 


æ Experiment 
Simulation 


— — 1D-Analysis 


0.8 


0.6 4 
0.4 4 
0.2 A 
0.0 1 1 1 1 
200 400 600 800 1000 
Mkr ol Fo, 


464 
æ Experiment 
Simulation 
08 — — 1D-Analysis J 
0.6 aan 4 
iam — 
n Ma 
s 
0.4 az 4 
Ea 
2 T =210°C 
4 all 
g 
02k 4 
0.0 L L 1 J 
200 400 600 800 1000 
m,kro lFioPs 
1.0 = = 
® Experiment 
Simulation 
= — 1D-Analysis 


0.8 - 


0.2 H 


—99n0 
i. a 290 C 


1 1 
400 600 


1 1 
800 1000 


mk /FioPs 


Fig. 7. Comparison of the numerical results for the conversion efficiency of methanol with the experimental and analytical results for three wall temperatures: (A) 


Tw =210°C, (B) Ty =220°C, and (C) Tw =230°. 


The numerical results for the conversion efficiency of 
methanol, along with the experimental and one-dimensional 
analytical results, are shown in Fig. 7 as a function of the 
dimensionless parameter m skg o/F1,ops for three wall tempera- 
tures. There is general agreement among all of the results; the 
numerical results are slightly lower than the one-dimensional 
results. 

An accurate determination of the CO produced from reaction 
(2) is important because CO deactivates the Pt catalyzed anode 
of fuel cells. The analytical solution for the mass fraction of CO 
produced in the catalyst can be obtained from [9] as follows: 
nico = 8 s 
mco = —McoLp x 10 

PgU 


(28) 


where co is the mass concentration of CO (ppm) leaving the 
reactor, ky is the rate constant of decomposition reaction defined 
by (1 — £) pskp, and Mco is the molecular weight of CO. Fig. 8 
shows the results for the mass fraction of CO leaving the reactor 
for three wall temperatures. The numerical and the analytical 
results are in fairly good agreement with one another but they 
differ significantly from the experimental data [13]. At the high- 
est wall temperature (230°C), the predictions are much greater 
than the experimental data. 

It is also noted that for the conditions studied, the numerical 
results yield little variation for the molar fraction of each species 
in the radial direction. In addition, for the constant wall temper- 
ature condition of this study, there is little variation of the gas 
temperature in both the radial and axial directions. 
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Fig. 8. Comparison of numerical results for the mass fraction of CO leaving the reactor with experimental and analytical results in three cases of wall temperature: 


(A) Tw =210°C, (B) Tw =220°C, and (C) Tw =230°C. 
4. Conclusions 


The conversion characteristics of the steam reforming of 
methanol in a microreactor tube packed with a commercial 
BASF F3-01(CuO/ZnO/A1203) catalyst have been investigated 
by numerical simulation, considering both the heat and mass 
transfer of the species. The numerical results for the conversion 
efficiency of methanol agree with a one-dimensional analyti- 
cal solution [9] and with experimental data [13] over the range 
of operating conditions. From the numerical results, a correla- 
tion for the conversion efficiency of methanol over the range 
of operating conditions is presented as a function of the oper- 
ating temperature and a dimensionless time parameter which 
represents the ratio of the characteristic time of the methanol 
flow to the time for the chemical reaction of methanol. For the 
mass fraction of CO leaving the reactor, the numerical and one- 
dimensional analytical results are in fairly good agreement, but 


they differ significantly from experimental data. At the highest 
wall temperature (230°C), the numerical and analytical results 
are much greater than the experimental data. The numerical 
results yield little variation for the molar fraction of each species 
in the radial direction. In addition, there is little variation of the 
gas temperature in both the radial and axial directions for the 
constant wall temperature condition of this study. 
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